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ABSTRACT 

We present the results of a long-term observation campaign of the extra- 
galactic wind- accreting black-hole X-ray binary LMC X-l, using the Propor- 
tional Counter Array on the Rossi X-Ray Timing Explorer (RXTE). The obser- 
vations show that LMC X-l’s accretion disk exhibits an anomalous temperature- 
luminosity relation. We use deep archival RXTE observations to show that large 
movements across the temperature-luminosity space occupied by the system can 
take place on time scales as short as half an hour. These changes cannot be 
adequately explained by perturbations that propagate from the outer disk on a 
viscous timescale. We propose instead that the apparent disk variations reflect 
rapid fluctuations within the Compton up-scattering coronal material, which oc- 
cults the inner parts of the disk. The expected relationship between the observed 
disk luminosity and apparent disk temperature derived from the variable occupa- 
tion model is quantitatively shown to be in good agreement with the observations. 
Two other observations support this picture: an inverse correlation between the 
flux in the power-law spectral component and the fitted inner disk temperature, 
and a near-constant total photon flux, suggesting that the inner disk is not ejected 
when a lower temperature is observed. 

Subject headings: Accretion, X-rays: Binaries (LMC X-l, LMC X-3, Cygnus X-l), 
Black Hole Physics, Stars: Winds 
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Introduction 


LMC X-l is (with Cyg X-l) one of only two known persistently luminous x-ray binaries 
consisting of a black hole accreting the wind of a massive blue star. The 10.91 ± 1.41M 0 


black hole is in a 3.9 0917 ±0.00005 day orbit flOrosz et al 


([Cowle y et al. 


20091 ) about an 07 III companion 


1995|). As the companion both drives a strong wine 


its Roche lobe, wind accretion feeds the black hole 


which accretes at an average of 0.1 6 L^^ flGou et ah 


low/hard state (Wilms et al 


2001 ]) 


Gou et al 


Nowak et al. 


and is far from filling 


2001;). The system, 


2009) has never been observed in the 


(120091 ) used x-ray data, including some of 


the RATA observations we use here, to derive a spin parameter a = 0.92(+0.05, —0.07) for 
the black hole in LMC X-l using fits to the disk blackbody component of the spectrum. 

LMC X-l’s persistent occupation of the high/soft state contrasts with the behavior 


of Cyg X-l , which it closely re sembles in other res pects. Cyg X 


black hole (IHerrero et ah 


companion (IBolton 


19951 ) in a 5.6 day orbit (IPaczvnskil 1 19741 ) about an 09.7 lab 


harbors a 10 M e 


19721 ). Though Cyg X-l is al so a w in d-accr eting system, its typical 


luminosity in the high state is ~ 0.03 Lem (Gierlinski et al. 


1 9991) . which is much lower than 


LMC X-l’s. This implies that the black hole in LMC X-l accretes at a higher rate than 
its counterpart in Cyg X-l, perhaps because the companions’ wind speeds and densities 
differ, or that accretion mechanisms with different radiative efficiencies are taking place. In 


l ow/ 


975 


-1 exhibits we. 

-documented transitions between the high 

i/soft state and t 

(Aerawal et al. 

1972; 

Baitv et al. 

1973; 

Heise et al. 

1975; 

Holt et al. 

1976 


Sanforc 


( Nowak et al 


et a. 


2001 


19751) . L MC X-l, in contrast, displays persistent soft-state emission 


Wilms et al. 


2001 ). While both LMC X-l and Cyg X-l’s softest states 


can be fitted with a power law plus disk blackbody model, they also differ markedly. The 


traditional soft state model, in whic 
accurately fits LMC X-l’s spectra I 


r the disk blackbody dominates the e nergy spectrum, 


Ebisawa et al. 


1989 


Yao et al. 


20051) . Cyg X-l’s 
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spectra, however, are energetically dominated by the power law in both the hard and soft 
states. Tims, while these two systems have similar companions and orbital properties, they 
differ in their spectral characteristics and evolution. 

An even more interesting comparison can be made between LMC X-l and LMC X-3, 
a black-hole binary that accretes via Roche-lobe overflow at comparable luminosity, as a 
means of highlighting the observational differences between disk and wind accretion. Much 
of this paper will be devoted to pointing out and interpreting these differences. 


LMC X-3 is a persistently bright x-ray binary in the Large Magellanic Cloud, and is 
the only other system with a dynamically confirmed black hole. As such, it provides a useful 


counterpoint to the win d-accre ting blac 


is_a 9 M 0 black hole (Ivan der Kli s et al . 


c hole binary systems mentioned above. LMC X-3 


19851 ) with an evolved B5 companion (ISoria et al. 


20011 ). A detailed listing of the system’s dynamical properties, along with those of LMC X-l 


and Cyg X-l, is provided in Table Q] Differences between the systems are expected to 
reflect differences in the accretion flows resulting from wind versus Roche-lobe overflow 
accretion: LMC X-3 accretes high angular momentum material through its first Lagrange 
point, while LMC X-l accretes low angular momentum material from its companion’s 
stellar wind. LMC X-3 and LMC X-l share spectral signatures of a persistent accretion 


disk, thoug 

i LM< 

Dunn et al. 

2010) 


C X-l’s disk spectrum defies the L oc T 4 relation (iMitsuda et al. 


1984; 


This paper examines how LMC X-l’s spectral properties and evolution fit with, 
and possibly extend, existing black hole binary accretion models. In particular, our 
results ultimately indicate that LMC X-l’s spectral behavior is consistent with sporadic 
obscuration of the innermost part of a stable accretion disk. 


Section 2 of this paper presents the data processing and spectral fitting techniques 
employed in obtaining the results presented in Section 3. Section 4 discusses how LMC X-l’s 
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Table 1:: System parameters for LMC X-l, LMC X-3, and Cygnus X-l, with references 
in parenthesis. The Eddington luminosity for each source was calculated from black hole 
mass. 



LMC X-l 

LMC X-3 

Cyg X-l 

Distance 

(kpc) 

48.10 ±2.22 (1) 

52.0 ±0.6 (2) 

2.5 (3) 

Black Hole 
Mass (M 0 ) 

10.91 ± 1.41 (1) 

11.1 ± 1.4 (4) 

10.1 (5) 

System 

Inclination 

(degrees) 

36.38 ± 1.92° (1) 

50.0° < i < 70.0° (12) 

36° < i < 67° (11) 

Companion 

Type 

07 III (8) 

B5 IV (6) 

09.7 lab (5) 

Companion 
Mass (M 0 ) 

31.79 ±3.48 (1) 

4.0 < M < 4.7 (6) 

17.8 (5) 

Companion 

Radius 

(* 0 ) 

17.0 ±0.8 (1) 

4.4 (6) 

22.7 ±2.3 (7) 

Orbital Pe- 
riod (days) 

3.90917 ±0.00005 (1) 

1.70479 ± 0.00004 (5) 

5.6 (3) 

Average 
Soft State 
Luminos- 
ity (% 

T 

~ 16% (8) 

< 30% (9) 

~ 3% (10) 

l^L'dZJ 





01 


( 2001 1 


Orosz et al.l (120091): 


(2009) 


(1198311 


; (5) 

; (9) 
; ( 12 ) 


Herrero et ah 


Nowak et al. 


21 Idi Benedettol (11997 1: 


1995); (6) 


Soria et al 


toiLO 


van derKlis et al 


11985) 


0 ) 


(3) 


Paczvnskil (11974 1 : (4) 


(1200111: (7) 


Gierlinski et ah 


Ziolkowski 


(1200511: (81 


(19991); (11) 


Gierl nski et al 


Gou et ah 


Davis fc Hartmann 
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unique spectral properties can be reconciled with existing accretion flow models, while 
Section 5 presents our conclusions and possible avenues for further investigation. 


2. Observations and Data Analysis 

The data analyzed in this work fall into two categories. The bulk of the observations 
come from our twice-weekly RXTE monitoring campaign of LMC X-l and LMC X-3. The 
LMC X-l data were complemented by long, individual RXTE observations available on 
the HEASARC. These data were reduced and analyzed using identical processes. The two 
classes of LMC X-l observations prove useful for probing different timescales of spectral 
variations. 


2.1. Observations 

The main observing campaign consisted of brief, uninterrupted observations of 
LMC X-l and LMC X-3, conducted twice a week since August 2007 and March 2006, 
respectively. These observations ranged in length between 506 and 5470 seconds, with 
average lengths of 1661 seconds for LMC X-l and 1877 seconds for LMC X-3. The 
LMC X-l observations were offset from the source by 15.764’ to reduce contamination from 
the nearby pulsar PSR B0540-69. These twice-weekly observations give reliable pictures 
of the typical X-ray spectra of LMC X-l and LMC X-3 over significant periods of time. 
While this is a major strength of these data sets, it also a weakness, in that they only have 
the capacity to reveal spectral variations that emerge over a minimum of several days. 

To examine shorter timescale variations in LMC X-l, we take advantage of the 
seventeen deep archival RXTE observations listed in Table [2j We “chopped” these data sets 
into shorter spectra (less than 90 min) corresponding to one orbit of the RXTE spacecraft. 
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These chopped archival observations can reveal spectral changes that emerge over fractions 
of a day. 


2.2. Data Analysis 


All data were reduced using standard HEASOFT v6.7 tools. Data taken within 
30 minutes after the satellite exited the South Atlantic Anomaly (SAA), or when the 
source was at less than 10° elevation above Earth’s limb, were excluded from subsequent 
analysis. All observations only considered data from the second proportional counter 
unit (PCU2) of RXTE’s Proportional Counter Array (PCA), due to its consistent 
performance throughout the mission. The mission-long faint background model 
pca_bkgd_cmf aintl7_eMv20051128 .mdl was used in this analysis, as the total count 
rate for both LMC sources lies below 40 photons/s/PCU. Data were extracted using the 
Standard2 data mode. 


We used uniform XSPE C v!2.5.1 (Arnaud 


Per RXTE recommendation ( Jahoda et al. 


1 9961 ) fitting procedures for all data. 


20061), systematic errors of 1% were added 


in quadrature to the Poisson counting noise errors. Because the detector response is 
poorly understood within the first three PCA channels, we exclude those channels from 
our analysis, resulting in a lower energy threshold just under 3 keV. Above 16 keV, the 
signal to noise ratio for both sources is too low to provide useful fits to the data, so our 


fits also excluc 


( Jahoda et al 


e_all data for energies over 16 keV. A 29% decrease in the PCA efficiency 


2006j) relative to an on-axis source resulted from the 15.764’ off-axis pointing 


of our LMC X-l observations. All fluxes derived from the off-axis LMC X-l observations 
were multiplied by a factor of 1.4 to account for this efficiency factor. Representative 


spectra, model fits, and fit residuals for LMC X-l and LMC X-3 are presented in Figures Q] 
and [21 respectively. 


Spectra were fitted with a disk blackbody plus power law emission model. The first 


spectral component is a multicolored blackbody spectrum generate d by t he accret ion c 


which is accounted for usin g the standard dis kbb XSPEC model flMitsuda et ah 


The simpl XSPEC model ( Steiner et al. 


198 41). 


isk, 


2009J) accounts for the power-law emission, which 


extends to high energies, and which is thought to be caus ed by Com pton up-scattering 


1976|). A pure 


of disk photons through a population of hotter electrons (jShapiro et al. 
power-law model runs the risk of promoting spectral fits that inaccurately trade disk 
luminosity for coronal energy flux. The XSPEC simpl model avoids this by turning over 
at low energies, and by convolving an arbitrary input photon spectrum with the specified 
Compton scattering prescription. It is, therefore, less prone to overestimating the amount 
of non-thermal emission. Simpl can account for either pure photon up-scattering, or for 
both up- and down-scattering. Because the average energy of coronal electrons is two orders 
of magnitude higher than the highest disk blackbody temperature, coronal emission from 
our sources can be safely approximated by pure up-scattering. 

Finally, to account for photoelectric absorption within the intervening ISM, we convolve 
the co mpound model with phabs. Our ap proa ch does not conflict with the findings of 


Levine fe Corbet 


(200Q) and 


Hanke et al. 


(120101) . which detect orbital phase variations in 
the absorption column towards LMC X-l. These orbital phase variations suggest that LMC 
X-l’s companion drives a significant wind. While the variable column density significantly 
effects LMC X-l’s observed soft x-ray emission, it does not impact the higher energy range 
examined in this analysis. 


Table [3] lists the fitting parameters and results from the following fitting procedure. 
The simpl power law index, the absorption column, and the inverse Compton scattering 


ratio normalized counts s _1 keV ^ 1 
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5 10 

Energy (keV) 


Fig. 1. — : normalFit and data-to-model ratio for LMC X-l data fit by the multi-colored 
disk and simpl (power-law) spectral model. 


ratio normalized counts s - 1 keV - 1 
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Fig. 2.—: 


Fit and data-to-model ratio for LMC X-3 data fitted by multi-colored disk and 


simpl (power-law) spectral model. 
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flag are all fr ozen a t the listed values. The absorption column densities are drawn from 


Nowak et al. 


(120011 1 for both LMC X-l and LMC X-3. The innermost disk temperature 
starts, but is not frozen, at the physically reasonable energy of 1 keV. The spectral fits 
output the inverse Compton photon scattering fraction, the innermost disk temperature, 
the overall normalization of the diskbb model component, and y 2 goodness-of-fit statistics. 
Once the model has been fitted to the data, it is used to calculate both comprehensive and 


component-specific photon and energy fluxes. 


The low photon counts from both LMC X-l and LMC X-3 limit the usable energy 
range in our 1.5 ksec observations to no higher than 16 keV. With counting statistics of this 
quality, it is not possible to simultaneously constrain both the power-law index and the 
inner-disk temperature. Since our observations cannot determine all six of the parameters in 
the phabs x simpl x diskbb model (see table), we freeze the power- law index, which is the 
most stable at long time scales. The deep archival observations of LMC X-l, prior to being 
chopped into orbit-by-orbit sections, have sufficient photons to simultaneously constrain all 
disk and power-law fit parameters. We use the average power law index derived from these 
long observations, which is consistent with a single value, to fix the power law index for 
all of the single-orbit LMC X-l spectra. For L MC X-3, we employ the average power law 


index found in the same way by 


Smith et al 


20071 ). In both cases, the values for all deep 


pointings are statistically consistent with the average. 


3. Results 

We find that LMC X-l remained in the soft state over the entire observation 
program, and has an anomalous disk-temperature-versus-luminosity relation relative to the 
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expectation of a disk blackbody. 


Figure [3] demonstrates how the luminosity of a typical accretion disk, such as the one 
present in LMC X-3’s soft st ate, follows the expected modified L oc Tf n Stefan-Boltzmann 


relation (IMitsuda et al. 


1984 1. The line in Figure |3l which represents this relation, was 


not fitted to the data beyond finding a suitable normalization constant. The relation 
nevertheless fits the data very well, aside from the highly uncertain points at low luminosities 
and high temperatures. These points come from observations taken after Feb. 23, 2009 
(MJD 54885), when LMC X-3 transitioned from a disk-dominated high/soft state to the 
low/hard state. The fitting routine we employ is not appropriate for the low/hard state, 
during which disk emission is suppressed and the power-law index hardens significantly. The 
poorly constrained values that appear in the lower right corner of Figure 21 well away from 
the Stefan-Boltzmann relation, result from using a high/soft state model to fit low/hard 
state data. 


Figure 0] similarly shows the disk temperature-luminosity relation for LMC X-l. The 
dashed lines are two normalizations of the same Stefan-Boltzmann relation used in Figure 
U and clearly cannot fit the data. 

Figure [3] compares how the inner disk temperatures of LMC X-3 and LMC X-l vary 
over time. The figure shows only a subset of the total observing campaign in order to better 
display short term variations from each source. The change in disk temperature from one 
observation to the next is markedly more continuous for LMC X-3 than for LMC X-l. 
Over longer time periods, however, LMC X-3’s inner disk temperature evolves significantly. 
Such long-term disk temperature evolution is absent in LMC X-l. 

Figure [3] shows how LMC X-l’s inner disk temperature and luminosity varied over the 
course of two of the archival observations listed in Table [2] The bottom panel shows data 
from the longest observation, while the top panel shows data from an observation with a 
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Inner Disk Temperature (keV) 


Fig. 3. — : The disk temperature-luminosity relationship observed in LMC X-3. The dashed 
line is the modified accretion disk Stefan-Boltzmann relationship. The points at low lumi- 
nosities and high temperatures, which lie far from the model, are due to a transition to the 
low/hard state towards the end of the observing campaign. 
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Inner Disk Temperature (keV) 


Fig. 4. — : The disk temperature-luminosity relationship observed in LMC X-l. The dashed 
lines are the modified accretion disk Stefan-Boltzmann relationship shown in Figure [3l and 
differ only in their overall normalization constants. 
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typical number of exposure intervals. 


The inner disk appears to evolve over a period of several hours. Rather than evolving 
along the Stefan-Boltzmann relation, the disk temperature and luminosity fall along a line 
of a very different slope. The order of the points is not uniform along this line either, 
suggesting motion faster than has been resolved by the data. 


Figure [7| shows LMC X-l’s variability in disk, coronal, and total luminosities, as 
well as in the total photon emission rate, over the course of our observing program. The 
mil d discrepancy be tween LMC X-l’s overall accretion luminosity as presented here and 


m 


Gouetah 


(120091 ) comes from the narrower energy range we employ in calculating the 
source’s overall flux. Both analyses, however, find similar disk-to-coronal flux ratios in LMC 
X-l. Percent variation, defined as standard deviation /mean, quantifies the extent to which 
the values in each panel vary. The system’s disk luminosity shows 14.5 % variation, while 
its coronal up-scattering luminosity shows 51.8 % variation. The system’s overall x-ray 
luminosity has 15.2 % variation. The total photon flux from LMC X-l is the system’s most 
stable feature, with only 12.5 % variation. While the total number of photons emanating 
from the system remains relatively steady, the amount of energy that the corona contributes 
to the total emission varies significantly. 


Figure [8] displays the relation between LMC X-l’s inner disk temperature and the 
inverse Compton scattering fraction. While observations with cooler disk temperatures and 
higher scattering fractions suffer larger fitting uncertainties, a Pearson correlation analysis 
confirms the existance of a statistically significant negative correlation between these two 
parameters, even when the two upper-leftmost points in Figure [8] are removed. Since the 
scattering fraction maps directly to the inverse Compton optical depth, Figure [8] shows that 
the inner edge of the disk appears cooler when there is more coronal material. 
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4. Discussion and Conclusions 


LMC X-l’s apparent defiance of the Stefan- Boltzmann relation may stem from 
variations in the disk itself, or from variations in the scattering medium surrounding it. 
Because the total number of photons is more nearly constant than the disk luminosity (see 
Figure [TJ), variations in the number of photons diverted to the power law by the corona 
seem more likely to explain the anomalous temperature/luminosity relations than variations 
within the disk itself. 


The first step in quantitatively distinguishing between these two explanations for the 
anomalous disk relation is to compare the observed disk variation timescales to timescales 
inherent to the disk and to the corona. In Figure [HI LMC X-l ’s disk randomly samples 
a wide swath of the temperature parameter space in the space of a few hours. Individual 
points in the archival observation are separated by an average of 2320 s, which provides an 
estimate of the timescale for apparent disk variations. We can compare this value to the 
disk’s viscous timescale to assess whether the disk’s structure is capable of changing over 
such a brief interval of time. 


The disk’s viscous timescale scales with both radius and viscosity, as (Frank et al 


2002 !) 


R 2 

tvisc ~ 

V 


( 1 ) 


The circularization radius, R C i rc , of the accreted wind material gives a lower limit for the 
outer disk radius. The disk radius, and therefore the viscous timescale, may actually be 
larger due to outward angular momentum transfer. Since we are primarily concerned with 
finding a lower limit on the viscous timescale, however, R c i rc provides a conservative value 
of the disk radius. The circularization radius is defined as the radius at which the accreted 
matter’s angular momentum due to Keplerian motion equals the an gular momentum it 


carried upon initial capture by the black hole. Following 


Frank et al. 


( 20021 ) we assume that 
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the wind’s initial angular momentum about the black hole goes as 

l ~ -a oR 2 acc (2) 

where R acc represents the black hole’s gravitational capture radius, and oj is the black hole’s 
orbital angular velocity about the companion. From there, we obtain 

G z Ml H u 2 


Rr.ire 


J rel 


( 3 ) 


where Mbh represents mass of the black hole, and 


2 2,2 
^ rel "h V or b • 


( 4 ) 


Therefore, the disk’s size is strongly dependent on the 


undergoes line-driven acceleration flCastor et al. 


1975 


speed o 


Owocki 


the stellar wind, wh ic 


1994- 


2000 1. The win d’s ve locity increases with distance from the star (ILamers & Casinelli 


Kudritzki fc Puls 


1999 


Kudritzki et al. 


1 989 ) as 


v w {r) = v a 


1 - 0.9983- 


R* 


( 5 ) 


Here, is the wind’s terminal velocity and is the companion’s radius. The factor of 
0.9983 is chosen such that v w equals the escape velocity at the star’s surface. The values 


of Voo and (3 are poorly constrained for O-star w inds, but 


km/s and from 0.5-1. 5, respectively (IKudritzki et al.ll 19891) , (ILamers fc Casinelli 


t ypically range from 


000-2000 


19991 ). We 


conservatively assume = 1700 km/s, in order to 


that consider clumping in O-star winds (|Ducci et af 


compare our results with other works 


2009)). We carry the effects of the full 


range of (3 values through all subsequent calculations to obtain a range of results. 

The final step in calculating R C i rc is determining v w where the wind reaches the 
accretion cylinder. To find R acc , we equate the wind’s gravitational and kinetic energy 
terms with respect to the black hole and find 

2 GMbh 


Rnrr 


-rel 


( 6 ) 
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Solving (J3J), (J5]) and dHJ) numerically over the full range of (3 values constrains v w to lie 
between 154 and 635 km/s by the time it reaches the black hole’s accretion cylinder. R c i rc 
is thus limi ted to values between 6.7 x 10 8 and 2.4 x 10 10 cm, in agreement with previous 


work (INowak et ah 


20011 ). 


Having found R acc , we turn to calculating the disk’s viscosity, v. The Shakura-Sunyaev 
ct-disk prescription posits 


v = ac.H 


where c s and H are the disk’s local sound speed and 


(7) 


of transient x-ray binary outbursts indicate a ~ 0.1 ([King et al. 


height, respectiy ely. Observations 


20071). Like LMC X-l, 


x-ray transients while in outburst host ionized accretion disks around a central black hole. 


These similarities motivate us to use a = 0 


1_ in the following calculations. Substituting 


equilibrium thin-disk formulae (IFrank et ah 
dZj) gives 

v = 1.8 x 10 14 a 4/5 iJ 10 m[ 1/4 fiJ 4 


20021 ) for the sound speed and disk height into 


Rmir, 

~TT 


1 / 2 ' 


3/10 


cur/g. 


( 8 ) 


Here, m\ is Mbh in units of M 0 , M±q is the accretion rate in units of 10 16 grams per second, 
Rio is radius in units of 10 10 cm, and R m i n is radius of the innermost stable orbit. Simple 
geometry requires that the accretion rate depends on the companion’s mass loss rate as 

2 


Mbh = - 


1 Rn 


4 \R 


^sys 


M * g/s 


(9) 


with a typical O-star M* of 10~ 5 M 0 per year. The factor of four in the denominator 
reflects the the fact that the black hole presents a circular cross-section to the spherically 
expanding wind. Equation (JHD shows that the accretion rate onto the black hole ranges 
between 2.2 x 10 1 ' and 9.6 x 10 18 g/s, depending on the chosen value of /3. These values are 
in good agreement with accretion rates derived from the observed x-ray luminosity using 


L ~ 0.1 c~Mbh-i 


( 10 ) 
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which range between 8.4 x 10 1 ' and 2 .0 x 10 18 g/s. The 0.1 efficiency factor in (fTOl) 


is recommended by 


Frank et al 


(120021 ). Equation [TUI arises from the condition that all 


gravitational potential energy up to the last stable orbit is radiated away. 

Using (HJ) to combine the results of (JHJ) and ([3]) gives 5 x 10 4 < t V i SC < 4 x 10 6 seconds. 
These values of t visc are at least an order of magnitude longer than the ~ 10 3 second 
timescale of inner disk variation indicated by Figure [6l This provides the first piece of 
evidence that the observed disk luminosity variations are not caused by a process inherent 
to the disk itself. But the calculation of the disk viscous timescale above assumes that 
the entire disk is involved in any change. The Roche-lobe accreting black-hole binary 
GRS 1915+105 shows very fast spec tral changes in many of its variability states near the 


Eddington luminosity (IBelloni et al. 


2Q00J), which appear to be repeated ejections of the 
innermost parts of the disk. But when the inner disk is ejected and the spectrum turns hard 
in GRS 1915+105, the total x-ray count rate also drops dramatically. It is the absence of 
any such change in count rate in LMC X-l that leads us to favor obscuration by marginally 
optically thick, ionized coronal material to an ejection mechanism. 


The slope of the temperature-luminosity relation that the chopped archival observations 
in Figure [6] occupy provides the second indication that LMC X-l’s apparent disk luminosity 
variations are caused by a process unrelated to the accretion disk. If the variations in disk 
luminosity were due to changes in the accretion rate, the points would move sequentially 
along the Stefan-Boltzmann relation shown in Figures [3] and [4] 


Since the anomalous disk temperature-luminosity relation is poorly explained by 
changes within the accretion disk, we ask whether changes in the corona could explain the 
relation. First, we consider which coronal properties are compatible with the observations. 
Then, we speculate how such a corona might form. 

Figure [8] provides a clue as to how the coronal material factors in LMC X-l’s anomalous 
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disk temperature-luminosity relation. The fraction of disk photons that suffer inverse 
Compton scattering while traversing the corona is lowest when the apparent inner disk 
temperature is around 0.9 keV. This value lies at the higher end of the observed range of 
disk temperatures for this source. The inverse Compton scattering fraction increases with 
decreasing inner disk temperature, up to a scattering fraction of nearly unity at Ti n of 
~ 0.4 keV. Such behavior is consistent with a steady disk whose inner radii are sometimes 
obscured by a cloud of energetic electrons. In that case, the number of disk photons, which 
form the seed photons for inverse Compton scattering in the overlying corona, would remain 
constant. The fraction of original photons that suffer up-scattering would depend only on 
the amount of coronal material present. Figure [7] shows that LMC X-l’s total photon flux 
remains largely constant, while its coronal luminosity varies substantially. 


To quanti fy t his model, we integrate the emission from the a— disk flux prescription 


(Frank et al. 


20021) outwards from inner radii at temperatures ranging between 0.4 and 


0.9 keV. Figure [9] shows the resulting temperature-luminosity relation, as fitted to all of 
the chopped archival datasets with four or more sub-intervals. The disk temperature and 
luminosity are constant in each observation. The amount of inner disk obscuration is the 
only variable that determines where the chopped observations fall along the best-fit line. 
The model appears consistent with all the data, particularly with observation 30087-01-04, 
which provides the strongest constraint. Note that these individual chopped observations 
run roughly perpendicular to the Stefan- Boltzmann relation in temperature-luminosity 
space. Suzaku or NuSTAR observations may allow the detection of variations of the 
power-law index on relevant (hour) timescales; hardness variations might then be interpreted 
as variations in Comptonizing optical depth, in which case hardness might be expected to 
correlate with the flux in the tail. 


Figures [9] and [8] raise questions about the analysis in 


Gou et al. 


020091 ). which derives 
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the spin of LMC X-l’s black hole from measurements of the inner disk temperature. 
Their analysis requires that the inner disk always be nnobscured. Onr analysis, however, 
shows that the degree of inner disk obscuration in LMC X-l varies significantly on short 
timescales. In systems such a s LMC X-3, where the disk is not obscured, there is no 


difficulty with the approach of 


Gou et al. 


( 20091 ) 


Why, in fact, does LMC X-3 have less apparent inner disk obscuration than LMC X-l 
does? We can imagine at least three possible answers. The first option is that the coronal 
optical depth is much larger, in general, in LMC X-l than in LMC X-3. This possibility is 
ruled out by the overlap between the two sources in terms of the flux ratio between the disk 
and power-law components; when the disk luminosity in LMC X-3is in decline, a significant 
power-law tail can form in that syste m without disrupting the excellent Stefan-Boltzmann 


relation for the disk component (ISmith et al 


2007|). A second possible answer supposes that 


the geometries of the coronae are very different: LMC X-3 has a geometrically large corona, 
while LMC X-l has one that is more centrally concentrated, so that it obscures more of the 
inner disk while producing the same amount of upscattering. A third possibility supposes 
that the corona in both cases has a conelike shape (possibly a jet, but not necessarily, from 
our evidence alone), and that the apparent differences between the two sourc es are due to 


an in clination effect. In LMC X-l, which has an inclination of 36.38 ± 1.92° (Orosz et al. 


c as we l ook dow n through it, 


19851) . the 


20091 ). the conical corona may block the inner parts of the disl 
while in LMC X-3, with an inclination between 50° and 70° (Ivan der Klis et al. 
cone presents itself in semi-profile, where it doesn’t obscure the inner disk but we can still 
see the Compton upscattered flux it produces coming out its sides. 


While either of the latter two pictures is consistent with our data, the last is more 
appealing for two reasons. Firstly, it does not require an unexplained difference in the 
form of the corona in the two binaries. Secondly, it may someday be possible to verify 
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or to disprove using x-ray polarization observations. A ’down the barrel’ jet observation 
(e.g. LMC X-l) would sport a lower over-all scattering polarization fraction than an 
observation of the same type of geometry viewed ’in profile’ (e.g. LMC X-3) would. The 
precise magnitude of the expected polarization fraction difference will govern how easily 
its presence or absence can be verified. In any case, future observations with the GEMS 
x-ray polarimeter will at least begin to explore the relevant paramet er space, prov ided the 


instrument can spend upwards of 10 5 seconds observing each source (jJahodal l2010h 


Finally, we turn to the question of how a population of energetic electrons forms and 
dissipates above the innermost region of the disk on timescales shorter than t V i SC . Figure 
[10] illustrates one possible accretion scenario that could account for LMC X-l’s unusual 
properties. 


In this picture, the incoming wind material sheds its net angular momentum in the 


“post shock” before accreting onto the black hole flHovle fe Lvttleton 


1939 ). Most of 


the shocked gas remains in the disk plane, and enters the accretion disk (middle arrow). 

A smaller portion of the heated post-shock matter launches on randomly distributed 
radial trajectories towards the black hole (upper and lower arrows). This component of 
the accretion flow rapidly heats on its radial descent into the black hole’s gravitational 
potential well. Unlike the disk, which is dense enough to cool efficiently via bremsstrahlung, 
the diffuse infalling matter cannot efficiently cool via free-free interactions. By the time 
the radially accreting matter passes over the inner disk, it can have reached electron 


'Nar avail & Yi 

1995; 

Ichimaru 

1977; 

Rees et ah 

1982; 

Esin et al. 

1998) 


19981 ). This model allows 


the corona to change independently of the disk. 


The presence of two accretion flow components, which are largely independent beyond 
the post-shock, is consistent with the constant observed total photon number and minimal 
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disk luminosity variation in LMC X-l. Changes in the wind mass flux, due perhaps to 
wind clumping, register immediately in the radial, coronal flow. Viscous diffusion in the 
disk, meanwhile, smooths out those same sharp mass discontinuities. While the coronal 
surface density, and therefore the total energy of coronal emission changes rapidly, the disk 
emission evolves more slowly and less dramatically. The scenario of independent flows with 
the thin disk remaining intact beneath the corona, even when the latter is optically thick, 
is consistent with the picture derived from observations of the hard state and hard-to-soft 


transitions in IE 1740.7-2942 and GRS 1758-258 ( Smith et af 


200211 . Recent observations of 


dativistic line profiles 

(Miller et al. 

2006 



2009). 


indicating that the inner disk remains present even in the hard state 
until extremely low luminosities, below 1% of Eddington, are reached (iTomsick et al. 

In the luminosity range of LMC X-l, we therefore expect the inner disk to remain present, 
whether the thin disk is actually continuous or whether its inner portion is recondensed 


from a disk that is disrupted at intermediate radii ([Liu et al 


201111 


Differing patterns of x-ray variability (e.g. Figure [5]) show promise as a means of 
identifying black hole companions when optical observations cannot. Fast changes in the 
inner disk temperature (Figure El), and an anomalous temperature/luminosity relation 
(Figure 0]), both combined with a stable net photon flux (Figure [7J, may signal the presence 
of a wind accretor, perhaps combined with a low inclination, just as long-term hysteresis in 
state c h anges appears to be a signature of Roche-lobe overflow accretion and a large disk 


(Smith et al. 


2002 


20071). 


The type of two-component accretion flow proposed in this paper is best suited to black 
hole binary systems where the companion drives a high-velocity stellar wind. At the same 
time, this particular model leaves ample room for other coronal generation mechanisms that 
have already been proposed for black hole binary systems accreting via Roche-lobe overflow. 
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At least two types of observations can either help to confirm or to disprove the 
hypothesis we have p roposed. The first focuses on recent x-ray observations of IC 10 X-l 


and NGC 300 X-l ( Barnard et al 


20081 ). Both systems have disk to power- law luminosity 


ratios, as well as high-mass main-sequence companions, that resemble LMC X-l’s. Binary 
systems that contain both an O-type star capable of driving such a wind and a black hole 
companion are rare, which has made this particular combination of accretion mechanisms 
particularly hard to observe. Further observations of IC 10 X-l and NGC 300 X-l, 
however, can potentially determine how deep their apparent similarities to LMC X-l 
run, and will provide more opportunities to test our model’s veracity. The second test 
relies on observations of intermediate state, low-disk-fraction black hole binaries. These 


systems ex hi 


(Dunn ct al. 


jit un usual temperature-luminosity relations akin to the ones seen in LMC X-l 


20111) . One can confirm whether the same mechanism is at work in both 
LMC X-l and an intermediate state system by measuring the variability (or lack thereof) 
of the latter’s total photon flux. Stable photon fluxes would support the proposition that 


these systems also contain stable, but variably-occulted, inner disks. 


LMC X-l has provided specific evidence that stellar mass black hole binaries may have 
x-ray spectral features which are unique to wind accretion. In a more general sense, the fact 
that LMC X-l shares certain x-ray spectral properties with other types of black hole binary 
systems motivates further investigation of its seemingly unorthodox accretion mechanism’s 
full scope and utility. 


This work was supported by NASA grant NNX09AC86G. The authors thank the 
anonymous referree for his or her constructive and insightful suggestions. L.R. thanks E. 
Ramirez-Ruiz, R. Strickler, and J. Naiman for productive discussions. 
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Observation 

UT (yyyy- 

Start Date 

Total Expo- 

Sub- 

ID 

mm-dd) 

(MJD) 

sure (s) 

Observations 

20188-01-02-00 

1996-12-30 

50447.438 

9919 

4 

20188-01-03-00 

1997-01-18 

50466.347 

9941 

3 

20188-01-05-00 

1997-03-09 

50516.354 

10189 

4 

20188-01-06-00 

1997-03-21 

50528.004 

10625 

4 

20188-01-07-00 

1997-04-16 

50554.193 

11484 

3 

20188-01-14-00 

1997-09-09 

50700.695 

10174 

4 

20188-01-18-00 

1997-10-10 

50731.611 

11436 

3 

30087-01-02-00 

1998-01-25 

50838.667 

10017 

3 

30087-01-03-00 

1998-02-20 

50864.370 

9948 

4 

30087-01-04-00 

1998-03-12 

50884.495 

10007 

7 

30087-01-06-00 

1998-05-06 

50939.149 

11396 

4 

30087-01-07-00 

1998-05-28 

50961.101 

5966 

2 

30087-01-08-00 

1998-06-28 

50992.097 

9921 

4 

30087-01-09-00 

1998-07-19 

51013.909 

10094 

4 

80118-01-06-02 

2004-01-11 

53015.340 

10908 

3 

80118-01-07-00 

2004-01-12 

53016.258 

10788 

3 

80118-01-08-00 

2004-01-10 

53014.236 

6608 

2 


Table 2:: Summary of the deep archival RXTE observations used to constrain the power 
law index, and to find upper limits on the apparent timescale of disk luminosity evolution. 
The former used only the parts of the spectra above 12 keV, while the latter relied only on 
regions of the spectra below 12 keV. 
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LMC X-l 

LMC X-3 

Inputs (frozen) 

nH (cm 2 ) 

7.2 x 10 21 (1) 

3.2 x 10 20 (1) 

r 

2.68 

2.34 (2) 

Scattering Flag 

1 

1 

Outputs (average values) 

kT e (keV) 

0.90 ±0.06 

1.07 ±0.18 

fsc 

0.13 ±0.07 

0.08 ±0.17 

X 2 /d.o.f. 

0.79 ±0.24 

0.79 ±0.77 


Table 3:: Input parameters and results for the phabs x (simpl x diskbb) fits to LMC X- 
1 and LMC X-3 data, with citations in parenthesis. Uncertainties in the outputs correspond 


to single standard deviations of the values about their average. Fixing simpl’s scatterin g 


fla g at a va l ue o f 1 forces the fits to account for only up-scattering. (1) 


Nowak et al. 


( 2 ) 


Smith et al 


(J2Q01|); 


1 20071 ) 


> LMC X-3 



<L 5.430x1 0 4 5.447 x 1 0 4 5.463x10 4 5.480x10 4 

Observation Date (MJD) 


Fig. 5. — : The innermost disk temperature plotted as a function of observation date for 
LMC X-3 (top) and LMC X-l (bottom). 
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Fig. 6. — : LMC X-l inner disk temperature vs. luminosity over the course of two repre- 
sentative archival observations. The top panel shows observation 20188-01-05-00, while the 
lower panel shows observation 30087-01-04-00. Black points, representing subsections of the 
archival observation, are numbered in chronological order. Gray points represent the data 
from our own observing campaign, and provide a basis for comparison. 
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Disk Luminosity vs. Date 



5.44x 1 0 4 5.46x1 0 4 5.48xl0 4 5.50x10 4 


Date (MJD) 

Fig. 7. — : From top to bottom: disk luminosity, coronal luminosity, total photon emission 


rate, and total luminosity of LMC X-lvs. observation date. Luminosities are shown in terms 
of the source’s Eddington luminosity. 
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Fig. 8. — : The fraction of total disk seed photons that undergo Compton np-scattering 
versus the inner disk temperature. 
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Fig. 9. — : The obscured inner disk model, fitted to all of the chopped archival observations 
containing more than three sub-intervals. Observation numbers are listed in the upper left 
corner of each plot. Black points represent the chopped observations, while gray points 
represent our data. The solid black line plots the obscured inner disk model, which has been 
normalized to fit the chopped observations. 
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Fig. 10. — : A cartoon of the accretion geometry of wind accretion systems. Incoming mate- 
rial from the companion’s wind (dashed gray lines) must dissipate its angular momentum in 
a shock front behind the compact object, a.k.a. the ’post shock’, before it can be accreted 
(thick arrow). 
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